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ABSTRACT 


The  effects  of  turbulence  on  the  performance  of  imagers  or 
on  beam  forming  optical  systems  are  well  expressed  by  the 
optical  transfer  function  (OTF)  or  its  magnitude,  the 
modulation  transfer  function  (MTF).  It  has  been  shown  that  the 
MTF  can  be  expressed  in  terms  of  the  Fried  model  by  means  of  a 
single  number,  the  turbulence  structure  constants  for  optical 
index,  C^,  provided  that  a  properly  path  weighted  value  is 
obtained.  Based  on  current  theory  the  weighting  needed  in  the 
path-position  weighted  value  of  is  a  function  of  (z/z0)5^3, 
heavily  emphasizing  the  region  near  the  imaging  system  or  beam 
forming  optics.  It  is  the  conclusion  of  this  report,  that  the 
path  weighting  function  (Z/Z0)5/3  is  correct. 
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I .  INTRODUCTION 


The  quantity  C^,  which  varies  with  the  micrometeorology, 
can  be  determined  in  its  properly  path-weighted  form  with  the 
slit  scanning  telescope  system  developed  at  NPS  by  Crittenden 
et  al  (Crittenden,  1978).  C^  can  also  be  computed  from 
measurements  of  the  temperature  structure  constant,  C*..  Point 
measurements  of  C^  are  based  on  assumptions  of  the  statistical 
form  of  the  turbulence,  while  optical  determination  of  c£  with 
the  slit-scanner  is  independent  of  this  statistical  nature. 

Anomolies  in  the  distribution  of  C^  due  to  air  flow  around 
obstructions  near  the  optics  is  a  common  problem.  This  has  a 
direct  naval  application  since  for  ship  mounted  optics,  the 
heavily  weighted  part  of  the  path  is  the  shipboard  end. 

The  objective  of  this  work  is  to  design  and  conduct  an 
experiment  to  determine  the  path-weighted  values  of  Czn  and 
compare  these  results  to  the  existing  model  for  the  path 
weighting  of  the  MTF.  The  model  used  is  that  proposed  by  the 
theoretical  treatment  of  Fried  (Fried,  1966)  which  predicts  a 
(Z/Zq ) 5 ^  dependence  where  ZQ  *  range. 


II .  THEORETICAL  CONSIDERATIONS 


A.  GENERAL 

Turbulence  and  temperature  gradients  in  the  atmosphere 
result  in  random  fluctuations  in  temperature  about  an  average 
value.  Similarly,  turbulence  and  humidity  gradients  result  in 
humidity  fluctuations.  The  optical  index  of  refraction  depends 
on  both  the  temperature  and  the  humidity,  so  that  the  index,  n, 
undergoes  fluctuations  around  its  mean  value. 

A  quantity  useful  in  characterizing  optical  turbulence 
effects  is  the  average  of  the  square  of  the  difference  of  the 
index  of  refraction  at  two  points  r  and  r  separated  by  a 
distance  r ,  i . e. , : 

<  [n(rt )  -  n(r2  )  ]2> 

For  isotropic  turbulence,  the  quantity  C^,  expressed  by 
Tatar ski  (Tatar ski,  1961): 

<[n(r,)  -  n(rz)]2>  =  r2/3  C2,  lQ  <<  r  <<  L0 

is  called  the  "index  of  refraction  turbulence  structure 
constant".  This  relationship  holds  within  the  "inertial 
sub- range", 

lQ  «  r  «  Lq, 

where  L  and  are  the  inner  and  outer  scales  of  turbulence 
o  o 

respectively.  These  scales  correspond  roughly  to  the  size  of 
the  smallest  eddies  of  turbulence,  about  a  millimeter,  and  the 
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largest  eddies,  often  taken  to  be  the  order  of  the  height  above 
the  surface  of  the  point  or  propagation  path  of  interest. 

The  quantity  is  a  property  at  a  point  in  the  atmosphere. 
The  mean  of  the  integrated  value  of  c£  over  the  optical  path, 
properly  path  weighted  becomes  the  primary  parameter  of 
interest  for  the  optical  properties  of  the  atmosphere  due  to 
turbulence . 

Values  of  Cn  at  specific  points  in  the  atmosphere  can  be 
measured  from  the  temperature  and  humidity  fluctuations.  The 
index  of  refraction  n,  is  related  to  the  absolute  temperature, 
T,  by 

n  -  1  *  79  x  10“®  p/T 

with  the  pressure  p  in  millibars.  The  average  of  the  square  of 
the  temperature  fluctuation  can  be  shown  to  follow  an  equation 
similar  to  that  for  index  fluctuation,  leading  to  a  definition 
of  the  temperature  turbulence  structure  constant  in  terms 
of: 


[T (f  )  -  T (?  )  ]*>  «  Ci  r2/3 
*  2  *- 

Using  the  expression  for  index  in  terms  of  p  and  T  above,  this 
leads  to  (Ochs,  et  al,  1969): 

Cft  -  (79  x  10'*  (p/T*)]2  C£ 
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Earlier  work  at  NPS  {Crittenden,  1978)  indicates  that  the 
contribution  to  the  index  variation  ^rom  humidity  fluctuations 
is  small  compared  to  temperature  fluctuations  and  is  therefore 
neglected  in  this  report. 

B.  TURBULENCE  EFFECTS  ON  OTF  AND  WANDER 

Turbulence  in  the  atmosphere,  as  expressed  by  a 
path- integrated  c£,  determines  the  performance  of  imaging  or 
beam- forming  optical  systems,  and  this  performance  can  be  used 
to  measure  the  path-integrated  C^. 

The  effects  of  atmospheric  turbulence  on  the  performance  of 
an  imaging  or  beam  projection  system  are  usually  treated 
(Fried,  1966;  Lutomirski  and  Yura,  1974)  by  considering  the 
atmosphere-plus-optics  combination  as  a  linear  system,  with  the 
image  being  the  response  of  the  system  to  the  object.  If  there 
were  no  sources  of  image  degradation  between  the  object  and  the 
image,  a  point  object  would  have  a  point  as  its  image. 

However,  diffraction,  caused  by  the  finite  aperture  of  the 
system,  is  always  present.  Within  the  present  context, 
refraction  and  diffraction  caused  by  the  index  of  refraction 
inhomogeneities  in  the  atmosphere  are  also  present.  These 
effects  cause  the  ideal  point  response  to  a  point  object  to  be 
smeared  out.  If  the  point  spread  function  has  the  same  shape 
regardless  of  its  position  in  the  image  plane,  then  the  image 
function  becomes  the  convolution  of  the  object  function  and  the 
point  spread  function.  In  this  case  the  convolution  theorem  of 
Fourier  transform  theory  yields: 
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i(vx,vy)  =  H(vx,vy)  •  o(vx,vy) 


where  ^vx»vy)»  H^vx' vy^'  and  °^vx'vy^  are  fc^e  two“^^mensional 

Fourier  transforms  of  the  image  function,  point  spread 

function,  and  object  function  respectively  and  the  values  are 

the  spatial  frequencies.  The  quantity  H(\>  ,  u  )  is  known  as  the 

x  y 

optical  transfer  function  (OTF)  and  its  magnitude  is  called  the 
modulation  transfer  function  (MTF).  The  optical  transfer 
function  is  thus  the  fundamental  quantity  which  cor.tains  the 
effects  of  the  atmopsheric  turbulence.  The  modulation  transfer 
function  is  also  used  to  characterize  the  effects  of 
turbulence,  but  it  lacks  information  on  phase. 

For  isotropic  turbulence  and  circularly  symmetric  optics, 
the  OTF  is  a  function  of  vxs  vy=  J  v)  .  In  this  case,  a 
one-dimensional  treatment  is  adequate.  Physically,  a 
one-dimensional  image  signal  is  obtained  by  scanning  the  image 
of  a  point  object  with  a  slit.  Then  the  point  spread  function 
is  replaced  by  the  line-spread  function: 


h(  x) 


**  Jh(x, 


y)  dy  and  i(v  )  =  H  fv  )  o(v) 


Theoretical  studies  (Fried,  1966;  Lutomirski  and  Yura, 
1971)  for  an  imaging  system  which  include  the  effects  of 
atmospheric  turbulence  show  that  the  total  optical  transfer 
function  can  be  expressed  as  the  product  of  the  transfer 
function  for  the  optics  and  the  transfer  function  for  the 


atmosphere. 
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Htot  M 


Hopt<v)  ’  Hatm<vl 


Because  of  the  stochastic  nature  of  turbulence,  quantities 
which  relate  to  it  can  be  understood  only  by  means  of  an 
ensemble  average.  The  problem  of  calculating  the  atmospheric 
MTP  is  considered  in  two  limiting  situations:  the  short  term 
or  short  exposure  (S.E.);  and  the  long  term  or  long  exposure 
(L.E.).  The  short  term  MTF  describes  images  obtained  in  the 
limit  of  time  intervals  sufficiently  short  that  the  turbulence 
can  be  considered  frozen.  The  average  short  term  MTF  is  then 
the  average  of  the  MTF 1 s  corresponding  to  many  such  images . 

The  long  term  MTF  results  from  taking  an  image  of  sufficiently 
long  term  which  sees  effectively  all  possible  turbulence 
configurations.  The  average  long  term  MTF  is,  therefore,  the 
same  as  that  from  a  single  long  time  integrated  exposure,  since 
every  turbulence  configuration  has  been  taken  into  account 
(properly  weighted). 

If  one  observes  a  series  of  line-spread  images,  one  would 
find  that  each  image  is  broadened  over  the  diffraction  limit 
imposed  by  the  optics  and  that  the  center  of  the  area  of  each 
image  appears  to  wander.  The  first  effect  is  called  image 
spread  and  the  second  image  wander.  The  short  exposure  MTF  is 
obtained  by  taking  the  average  (and  then  the  Fourier  transform) 
of  successive  line-spread  functions  in  which  each  line-spread 
function  is  shifted  so  that  the  center  of  the  area  of  each 
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curve  has  a  common  origin.  This  procedure  removes  the  image 
wander.  The  long  exposure  MTP  is  obtained  by  not  performing 
this  shift  in  origin,  therefore  including  image  wander.  The 
long  exposure  MTF  is  seen  to  be  the  result  of  image  spread  (the 
short  term  MTF)  plus  image  wander. 

Several  theoretical  models  have  been  proposed  from  which 
the  L.E.  MTF  can  be  derived.  The  only  theoretical  treatment  of 
the  S.E.  MTF  thus  far  is  that  of  Fried  (Fried,  1966).  The 
result  can  be  summarized  by  the  following  expression: 

M  (f )  =  <H(f)>  *  exp  £-57.64  (5  C2  zQ  X-1/3  f5/s  [1  -  a  <^) 1  /3]  j 

where 

X  »  wavelength 

f  *  angular  spatial  frequency  ( cycles/radian) ,  *  fv 
v  »  linear  spatial  frequency  {cycles/meter ) 

F  *  focal  length  of  system 
D  *  diameter  of  optics 
zQ  a  range 

10  for  the  L.E.  OTF 

1/2  for  the  far-field  S.E.  MTF,  D  <<  (zX)l/z 
1  for  the  near-field  S.E.  MTF,  D  >  (zX)1/2 

(1  for  a  plane  wave 
3/8  for  a  spherical  wave 


Use  of  this  quantity,  f,  the  angular  spatial  frequency 
appearing  in  this  expression,  avoids  the  necessity  to  know  the 


focal  length  or  magnification  of  a  system.  Most  experimental 
situations  which  employ  a  diverging  laser  beam  as  the  light 
source  or  object  can  be  considered  to  correspond  to  the 
spherical  wave  model.  Plane  waves  are  rarely  encountered  and 
will  not  be  considered  in  this  report. 

C.  PATH  POSITION  WEIGHTING  OF  c£  FOR  THE  MTF 

The  expression  (in  the  preceding  section)  for  the  MTF 

2 

applies  to  the  case  where  Cn  is  independent  of  the  position 
along  the  propagation  path.  If  Cn  varies  along  the  path,  thi 
equation  generalizes  for  the  L.E.  spherical-wave  MTF  to 

-zo  i 

M (f )  =  exp  £-57.64  X"1/3  f S/*J  C*  •  (z/z0)5/3dzj 

That  is,  is  replaced  by  a  weighted  mean  given  by 

Cn  *  <z/zo)5/3  dz 

The  experimentally  determined  MTF  can  be  used  to  infer  this 
weighted  mean  for  C^* 

The  above  equation  expresses  the  path  weighting  of  C* 
indicated  in  Figure  1-1  in  which  the  path  region  nearest  the 
imaging,  or  beam  projecting  system  at  z  =*  zQ  has  the  largest 
contribution.  The  opposite  end  of  the  path  has  zero 
contribution. 
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Detector 


III.  GENERAL  PLAN  OF  EXPERIMENTAL  PROGRAM 


A.  GENERAL  CONSIDERATIONS 

The  general  objective  is  to  confirm  the  path-weighting 
model  for  the  MTP.  If  the  background  in  a  controlled 
environment  is  several  orders  of  magnitude  below  the  of  a 
highly  localized  temperature  disturbance,  then  the  optically 
measured  will  be  due  primarily  to  the  injected  disturbance 
and  the  contribution  of  the  remainder  of  the  path  may  be 
ignored. 

The  major  requirements  of  the  experiment  are  then  to 
construct  an  environment  with  sufficiently  low  and  uniform 
background  and  to  construct  a  heat  source  which  adequately 
meets  the  isotropic  turbulence  requirements.  This  heat  source 
will  require  a  method  of  monitoring  to  ensure  a  nearly  uniform 
disturbance  as  it  is  repositioned  along  the  path  in  the 
controlled  environment.  If  slight  non-uniformities  exist,  they 
may  be  used  as  normalizing  factors  for  the  optically  determined 
values  of  C^. 

B .  PROCEDURE 

A  sixty-one  meter  tunnel  was  constructed  for  the  controlled 
environment  and  a  suitable  heat  source  designed  for  the 
temperature  disturbance.  By  placing  the  heat  source  at 
progressive  locations  in  the  "tunnel",  measuring  the  optically 
determined  values  of  C^  from  the  MTP  as  determined  by  the 
slit-scanner  and  on-line  data  processing  system,  the  results 
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were  plotted  as  versus  the  temperature  disturbance  location 
in  the  tunnel.  These  are  compared  to  the  theoretical  curves 
fitted  to  the  data  taken  on  the  same  graphs. 

Monitoring  of  the  heat  source  was  achieved  by  using 
platinum  temperature  probes  across  an  AC  bridge  to  measure  c£ 
and  thereby  compute  c£  at  the  location  of  the  temperature 
disturbance  location.  These  values  will  not  necessarily 
coincide  with  the  optically  determined  values  of  due  to 
their  dependence  on  the  statistical  nature  of  the  temperature 
of  fluctuations  of  the  heat  source  (micrometeorology) .  The 
optical  measurements  assume  a  uniform  over  the  path.  As 
mentioned  earlier,  the  controlled  environment  ambient  c£  -is 
much  less  than  the  cf.  of  the  heat  source. 
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IV.  MEASUREMENT  TECHNIQUES 


A .  MTF 

1 .  Long  Term,  Short  Term,  Wander 

The  performance  of  an  imaging  system,  or  of  a  beam 
projection  system  such  as  the  HEL  and  pointer  (tracker)  can  be 
characterized  and  predicted  in  terms  of  the  MTF  of  the 
atmosphere.  The  experimental  program  has  been  designed  to 
provide  verification  of  the  theoretical  model  for  the  path¬ 
weighting  of  the  MTF. 

As  indicated  in  Section  II,  the  MTF  is  expected  to  depend 
on  the  level  of  turbulence  as  expressed  by  the  path-weighted 
value  of  C^,  and  on  the  wavelength  and  range.  Two  experiments 
were  carried  out  to  determine  the  effect  of  a  temperature 
disturbance  on  the  path-weighting  of  the  MTF.  These  have  been 
carried  out  in  a  controlled  environment  as  free  from  other 
turbulence  as  practicable  with  uniform  values  of  Cn  over  the 
path  except  at  the  points  of  interest. 

Measurements  of  the  MTF  yield  long  term,  short  term,  and 
wander.  The  internal  consistency  of  these  measurements  are  to 
be  compared. 

To  calibrate  the  MTF  in  angular  spatial  frequency  with  a 
maximum  scale  of  500  microradians  consisting  of  1024  data 
points,  a  calibration  signal  is  injected  by  placing  a 
diffraction  grating  in  front  of  the  optics.  The  spacing 
between  first  order  peaks  is  known  to  be  142.6  microradians,  or 
A/d  »  142.6/2  »  71.3  microradians.  The  dwell  time  is  the 
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number  of  microseconds  per  point.  For  a  dwell  time  of  20,  the 
time  between  the  first  order  interference  peaks 

»  (20  x  10”6  usec/point)  (1024  points)  (142.6/500) 

»  5.841  milliseconds 

Using  a  Hewlett-Packard  1743  oscilloscope,  this  time  is  set  on 
the  calibration  signal  by  adjusting  the  gain  of  the  scanning 
mirror  voltage.  Center  display  of  the  signal  is  obtained  by 
adjusting  the  scanning  mirror  D.C.  offset. 

The  calibration  grating  is  then  removed  and  the  desired 
signal  observed. 

A  Wavetek  model  145  function  generator  provides  the  sweep 
signal  to  the  scanning  mirror  controller  and  a  TTL  output  to  an 
I EC  P-12  pulse  generator  which  provides  a  delayed  trigger  to 
the  data  processor  NIC-80  described  in  section  IV-2.  This 
delayed  trigger  prevents  the  NIC-80  from  seeing  the  sweep-back 
signal  from  the  scanning  mirror.  Figure  IV-1  diagrams  the 
optical  data  processing  equipment. 

2.  On-Line  Data  Reduction 

The  on-line  data  reduction  unit  is  a  Nicolet  Instrument 
Company  NIC-80  data  processor.  This  is  a  data  averager 
interfaced  directly  to  Hewlett-Packard  9825B  programmable 
calculator.  The  HP-982 5B  acts  as  the  controller  and 
input-output  device  to  the  NIC-80  and  controls  the  plotter  and 
printer. 

The  NIC-80  digitizes  the  analog  signal  from  the 
detector,  processes  the  raw  data  to  obtain  the  experimental 
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MTF,  wander  and  C^.  The  processed  experimental  data  is  then 
transferred  to  the  HP-982 5B  which  is  better  suited  for  floating 
point  calculations.  The  HP-9825B  calculates  the  best  fit  curve 
and  plots  out  the  results  in  final  form. 

3.  Optical  Equipment 

The  optical  equipment  consists  of  an  8  inch  diameter 
Celestron,  Schmidt  telescope  with  an  equivalent  focal  length  of 
160  inches  and  an  overall  length  of  about  50  centimeters.  It 
has  a  .344  central  obscuration  ratio.  Mated  to  the  image  end 
of  the  telescope  is  the  scanning  mirror  assembly.  It  consits 
of  an  optically  flat  mirror  mounted  on  a  servo  motor  to  provide 
the  sawtooth  scan  of  the  image  across  an  adjustable  slit.  The 
mirror  turns  the  optical  path  90  degrees.  The  scan  produces 
the  line  spread  function  as  a  function  of  time.  A  trigger 
signal  associated  with  initiation  of  each  scan  indexes  the 
position  of  the  image. 

The  detector,  mounted  after  the  slit  in  the 
slit-scanner,  is  an  RCA  silicon  photodiode  type  C30872.  It  is 
operated  at  370  volts  D.C.  with  a  9  volt  battery  powering  the 
source  follower.  The  detector  used  does  not  require  cyrogenic 
cooling  which  is  a  large  improvement  in  reliability  and  ease  of 
use  over  a  detector  that  did  during  earlier  trials. 

The  telescope  optics  are  supported  on  a  rigid  tripod 
base,  weighted  with  50  pounds  of  lead  which  prevents  mechanical 
vibration  of  the  telescope  system.  The  optics  and  detector  are 
sealed  as  a  part  of  the  tunnel  during  optical  measurements. 
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The  signal  from  the  detector  is  preamplified  by  a 
high-gain,  low-noise  PAR  113  amplifier.  The  signal  is 
processed  immediately  with  the  on-line  equipment. 

The  source  for  the  MTP  system  is  a  Helium-Neon  laser 
transmitted  from  the  far  end  of  the  tunnel.  The  laser  beam 
passes  through  two  polarizers  to  reduce  its  intensity  and  is 
diverged  by  a  10  centimeter  lens  to  provide  a  spot  size  at  the 
MTF  telescope  of  about  .5  meters  wide.  This  provides  uniform 
illumination  over  the  aperture  of  the  telescope  and  closely 
approximates  light  diverging  as  a  spherical  wave  from  a  point 
source.  The  MTP  telescope  then  images  this  point  source.  The 
laser  is  operated  in  the  CW  mode  for  this  purpose.  The  laser 
wavelength  utilized  for  MTF  measurements  is  the  red  HeNe  at 
0.6328  micrometers. 

4»  Controlled  Environment  (Tunnel) 

The  tunnel  was  constructed  of  three-eighths  inch 
construction  plywood  with  completed  sections  measuring  2  feet 
by  2  feet  by  8  feet  (61  cm  by  61  cm  by  2.44  m) .  25  tunnel 
sections  were  made  for  a  tunnel  length  of  200  feet  or  61 
meters.  Prior  to  construction  an  underseal  of  water  repellant 
was  applied  both  inside  and  out.  The  inside  was  then  painted 
with  flat  black  and  the  outside  with  an  exterior  white.  The 
individual  sections  were  constructed  with  12  metal  edge 
brackets,  8  end  corner  brackets,  and  56  wood  screws.  The 
joints  were  sealed  with  a  heavy  duty  tape  to  prevent  air  flow 
through  any  of  the  seams. 
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A  special  heat  source  tunnel  section  was  built  by 


modifying  an  existing  section.  An  opening  about  1  1/2  times 
the  area  of  the  heat  source  was  cut  in  the  top  at  one  end  and  a 
64  cm  high  chimney  constructed  over  it.  A  wooden  baffle  sheet 
was  placed  over  the  top  of  the  chimney  (supported  by  small 
blocks)  in  a  manner  to  allow  air  flow  but  prevent  down  drafts 
from  entering  the  tunnel.  A  wide  slit  was  cut  in  the  bottom  of 
the  tunnel  section  under  the  chimney  to  allow  for  air  flow 
across  the  heat  source  by  convection.  The  slit  was  baffled  to 
prevent  entry  of  side  drafts.  Finally  a  small  slit  was  cut  in 
the  side  of  the  tunnel  section  for  the  insertion  of  the 
platinum  temperature  probes  above  the  heat  source.  The  heat 
source  was  then  placed  under  the  chimney  in  the  bottom  of  the 
tunnel  section  (Figure  IV-2) . 

The  tunnel/  located  in  the  basembnt  of  Spanagel  Hall  at 
NPS,  was  assembled  for  the  experimental  runs  only  after  working 
hours  since  it  blocks  several  doorways  and  was  disassembled 
after  each  run.  When  assembled,  the  joining  and  sealing  of 
tunnel  sections  was  achieved  by  taping.  During  the  experi¬ 
mental  runs,  the  heat  source  was  injected  every  other  tunnel 
section  or  every  4.88  meters  along  the  61  meter  path. 

A  table  was  positioned  at  each  end  of  the  tunnel  to 
house  the  optics  and  detector  at  one  end  and  the  laser  and 
lenses  at  the  other.  During  the  runs,  these  tables  were  sealed 
and  became  an  integral  part  of  the  tunnel’s  controlled 
environment. 
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If  the  tunnel  were  a  perfect  controlled  environment 
with  no  temperature  disturbances,  a  point  source  at  one  end 
would  be  imaged  by  the  optics  as  a  point  source  at  the  other 
end  and  the  MTF  would  be  a  straight  line  (normalized  C2n  -  1 
independent  of  spatial  frequency) .  Figure  IV-3  shows  the  very 
high  degree  of  "perfection H  achieved  for  the  controlled 
environment  inside  the  tunnel. 

5.  Heat  Source 

As  discussed  in  Section  II,  a  region  of  isotropic 
turbulence  for  lQ  <<r<<  L0  is  required  (or  approximately  so) 
from  whatever  heat  source  would  be  used.  In  this  region  Cn 
must  be  reasonably  uniform.  For  several  heat  sources  the  C  as 
a  function  of  r  (the  probe  spacing  for  the  from  c£  equipment 
described  in  Sections  V-B  and  C)  was  plotted.  The  heat  section 
used  gave  the  best  results  and  the  data  is  displayed  in  Figure 
IV-4. 

The  heat  source  consists  of  a  two  conductor  heating 
element  bent  in  reversing  V's  mounted  on  a  wooden  stand  with 
ceramic  standoffs.  The  wide  spacings  of  the  elements  were  to 
assist  in  providing  the  larger  turbulence  size  which  was 
considered  to  be  about  the  height  of  the  probes  above  the  heat 
source,  22  cm  (Lo) .  A  screen  of  wire  spacing  .15  cm  was  placed 
over  the  heating  element  to  provide  the  smaller  turbulence  size 

U0)- 

Power  was  supplied  by  a  Variac  set  at  70  volts  AC  which 
produced  a  temperature  22  centimeters  above  the  heat  source  of 
about  27-30  Celsius. 
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Assuming  £0<<r0<<L0  or  .15  cm  <<r<<  22  cm  the  region 
of  isotropic  turbulence  would  be  expected  to  be  in  the  vicinity 
of  1.5  to  2.5  centimeters.  Figure  XV-4  shows  a  peak  around  2.0 
centimeters  the  remainder  of  the  curve  falling  off  as  r (the 
path  dependence  of  Cn) .  The  2  centimeter  region  was  considered 
to  be  approximately  isotropic.  A  probe  spacing  of  6 
centimeters  was  eventually  decided  upon  due  to  the  sharp  peak 

at  2  centimeters.  A  more  linear  section  of  the  curve  is 

2 

desired  for  more  uniform  results  since  the  measured  Cn  will  be 
used  to  monitor  the  stability  of  the  heat  source  and  provide 
weighting  for  the  optically  determined  results. 

Thus  decided  upon,  the  heat  source  was  placed  in  the 
bottom  of  the  special  tunnel  section  designed  for  it  (Figures 
IV  -  5  and  6) . 

B.  Cj.  POINT  DETERMINATION 
To  determine  C^, 

C|  =  <  [T  (r, )  -  T  (r2 )  ]  2>r“:^3 

T  must  be  measured.  The  value  of  r,  the  temperature  probe 
spacing  in  meters  has  been  picked  at  .06  meters  as  discussed  in 
the  previous  section. 

Using  two  platinum  wire  temperature  probes,  (  .0036°  K“  l) 

each  with  a  resistance  of  approximately  56  ohms,  across  an  AC 
Bridge,  the  output  is  calibrated  for  AT.  The  bridge 
sensitivity  was  set  and  measured  at  6.67  volts/ohm. 
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AT  per  volt 


(.0036  °K-1 ) (56  OHMS) (6.67  Volt/OHM) 


AT  per  volt  =  .744  K  per  volt. 

Using  a  spectrum  analyzer  on  the  AT  fluctuations,  the  half 
power  point  for  the  AT  fluctuations  was  determined  to  be  about 
8  hertz.  A  Hewlett-Packard  3438  A  System  Voltmeter,  sampling 
at  an  11  hertz  rate  for  500  samples  (45  seconds) ,  is  then  used 
to  determine  the  values  of  AT.  The  system  voltmeter  is 
controlled  by  a  Hewlett-Packard  9825B  programmable  calculator. 
It  squares  each  AT  sample  and  averages  the  500  squared  samples 
Multiplying  by  r-2^3  ,  the  calculator  outputs  c£. 

C.  C*  FROM  C2t 

,  2 

To  determine  Cn  , 

C2  [79  x  10~6  (P/T2) ]2  C| 

the  absolute  temperature  of  the  heat  source  must  be  measured. 

Using  a  Copper-Constantan  thermocouple  referenced  to  0° 
Celsius  (ice  water)  and  sampling  at  a  3  hertz  rate  using  a 
Hewlett-Packard  3490A  multimeter  controlled  by  the 
Hewlett-Packard  9825B  calculator,  forty  samples  of  temperature 
are  taken  and  averaged.  Twenty  are  before  the  sampling  for  Ct 
and  twenty  after. 

T  =*  < Thermo-Couple  Voltage>  32809.524  °C  per  volt.  Adding 


273.15  gives  the  temperature  in  degrees  Kelvin. 


Assuming  a  constant  atmospheric  pressure  of  1013.25 
millibars,  the  H/P  982 5B  calculates  to  be  used  for 
monitoring  the  heat  source  and  weighting  of  the  optically 
determined  values  of  C& 

Figure  IV-7  shows  a  block  diagram  of  the  equipment  used 
in  the  ,  T  and  measurements  and  calculations. 
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IV- 1.  Optical  Data  Processing  Equipment 


Figure  IV-4.  Isotropic  Turbulence  Region  Determination 


V.  RESULTS 


A.  EXPERIMENTAL  DATA 

1.  22  March  1983 

During  the  evening  of  22  March  1983,  the  first  experimental 
run  was  made  in  the  basement  of  Spanagel  Hall  at  NPS.  A 
"quiet"  tunnel  result  similar  to  Figure  IV-3  was  bbtained 
confirming  the  quality  of  the  tunnel's  environment.  The  point 
values  of  from  were  fairly  uniform  along  the  path  with 
an  apparent  systematic  increase  on  the  detector  end  of  the 
tunnel.  The  short  term  results  (Figure  V-l)  were  in  excellent 
agreement  with  theory  but  the  long  term  results  (Figure  V-2) 
though  good,  were  not  as  convincing.  The  results  from  wander 
(omitted)  were  totally  random  in  appearance. 

Since  wander  is  the  shift  in  image  center,  vibrations  could 
cause  these  erratic  results.  The  initial  telescope  mount  was 
not  the  one  described  in  Section  IV-A,  2.  The  originaL 
mounting  was  noted  as  being  susceptible  to  vibration  and  the 
wander  data  on  this  run  motivated  its  alteration. 

2.  11  April  1983 

After  the  telescope  vibration  problem  was  corrected,  a 
second  experimental  run  was  conducted  during  the  evening  of  11 
April  1983  in  the  same  location. 

from  point  values  (Figure  V-3)  were  very  uniform 
over  the  path.  Again  the  short  term  results  (Figure  V-4)  were 
in  excellent  agreement  with  theory  while  the  long  term  results 


(Figure  V-5)  were  not  as  good,  but  still  showed  the  general 
shape  of  the  curve.  The  wander  results  (Figure  V-6),  though 
vastly  improved,  cannot  be  used  to  substantiate  the  theoretical 
model. 

B.  ANALYSIS  OF  DATA 

The  following  conclusions  are  drawn  from  the  experimental 
data:  \ 

1.  The  short  term  results  confirm  the  Fried  model  for  the 
path-weighting  of  the  MTF. 

2.  The  long  term  results  give  good  credibility  to  the  Fried 
model,  but  cannot  be  used  to  positively  confirm  it. 

3.  The  wander  results  do  not  support  the  theoretical  model. 
The  overall  conclusions  is  that  the  Fried  model  is  correct 

and  that  systematic  noise  is  responsible  for  the  poor  quality 
of  the  wander  results  and  aberrations  in  the  long  term  results. 
As  described  earlier,  wander  is  not  removed  during  the 
processing  of  the  long  term  data  so  any  physical  or  electrical 
vibration  of  the  experimental  system  will  affect  the  wander 
primarily  and  the  long  term  to  a  lesser  degree. 

C .  RECOMMENDAT IONS 

Since  lead  bricks  weighted  all  physical  components  of  the 
system,  the  likelihood  of  mechanical  vibrations  during  the 
second  run  is  low. 

A  large  source  of  possible  error  seems  to  be  in  the  sweep 
generation  of  the  scanning  mirror.  The  flatness  of  the 
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sawtooth  sweep  voltage  was  adjusted  by  eye.  To  ensure 
linearity  over  the  sampling  region  of  the  sweep,  the  Fourier 
transform  of  the  ramp  function  could  be  computed  and  compared 
to  a  spectrum  analysis  of  the  scan-mirror  output  signal.  Any 
deviation  from  the  predicted  components  would  be  a  source  of 
mirror  vibration  and  hence,  wander. 
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